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Abstract 

In this paper we consider a reduction of a non-homogeneous linear system of first order operator equations to a totally 
reduced system, derived by using the characteristic polynomial of the system matrix. Variables in the new system 
are separated by higher order linear operator equations. 
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1. Introduction 



Let -RT be a field, V a vector space over K and \ei A -.V — > V he a, linear operator. Linear system of 
first order ^-operator equations with constant coefficients in unknowns xi, is: 



^(a;i) = 6iixi + 6122:2 + ■ . • + hinXn + ^\ 

A(x2)^h2\Xx + h22X2 + ■ . • + h2n,Xn + ^2 
^ A{Xn)^hn\Xx + hn2Xl + ■ ■ • + ferm^Jn + ^ 



(1) 



nxl 



for hij G K and tpi S V . We say that B = G iiT"^" is the system matrix, and <^= [<pi ... G V 

is the free column. If (^=0 system ([T]) is called homogeneous. Otherwise it is called non- homogeneous. 

Let x—\x\ ... Xr^ be a column of unknowns and A : — y ynxi j^g vector operator defined 
componentwise A{x) = [A{xi) . . . A{xn)]'^ ■ Then system ([1]) can be written in the following vector form: 

A{x) Bx + 0. (2) 

A solution of vector equation © is any column veV^^^^ which satisfies 
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Powers of operator A are defined as usual A* = A'^^o A assuming that A" = / : F — > V is tlie identity 
operator. By n*'* order linear A-operator equation with constant coefficients, in unknown x, we 

mean: 

A'^ix) + hA''-\x) + ... + hj{x) = if, (3) 

where 5i,. . .,6„ Gif are coefhcients and is non- homogeneous term. If <p = equation ([3]) is called homo- 
geneous. Otherwise it is non- homogeneous. A solution of equation ([3]) is any vector v&V which satisfies (|3]). 

2. Notation 

Let M be an n-square matrix. We denote by: 

0k) ~Sk{M) the sum of its principal minors of order k {1 < k < n), 

{d\.) 6\. =5jj.(Af) the sum of its principal minors of order k containing column {1 < i,k < n). 

Let vi, . . . ,Vn be elements of K. We write M^{vi, . . . , Vn) gX"^" for the matrix obtained by substituting 
v = [vi ... Vn]'^ in place of i"* column of M. Furthermore, it is convenient to use: 

6liM; v) = 61{M; v,,..., v„) = SUm^v^, . . . , v^)). (4) 

The characteristic polynomial Ab(A) of the matrix B<eK"'^^ has the form 

Ab(A) = det(A/ ~B)^ A" + rfiA"-i + . . . -I- d„_iA -I- d„, 

where dk = (-l)''5fe(S), 1 < fc < n; see [3, p. 78]. 

Denote by B{X) the adjoint matrix oi XI — B and let Bq, Bi, . . . , Bn-2, Bn-i be rt-square matrices over 
K determined by B{X) = adj(A/-B) = A^-^Sq + X'-^Bi + ... + AB„_2 + B^-i- 

Recall that (A/ - B)B{\) = Ab{\)I = B{\){XI-B). 

The recurrence Bo = I; Bk = Bk-iB + dkl for \ < k < n, follows from the equation 
B{\){\I - B) = 1^b{X)I, see [3, p.91]. 

3. Main Results 

The primary purpose of our paper is to reduce initial system ([2]) to a system of higher order linear 
^-operator equations, called totally reduced system. It has separated variables and symmetrical form. 
The following theorem will be useful in proving our main result. 

Theorem 3.1 Assume that the linear system oj first order A-operator equations is given in the form 
A{x) — Bx + ip, and that matrices Bq, . . . , -B„_i are coefficients of the matrix polynomial B{X) = adj(A/— _B). 
Then it holds: 

n 

(As(A))(x) = ^Bfc-i(l"-'=(^)). (5) 
fe=i 

Proof. Let Lb : 1/"^^ — > ynxi ^ Unear operator defined by Lb{x) — Bx. Replacing A/ by A in the 
equation Ab(A)/ — B{\){\I — B) we obtain 

As(i*) = o (1- Lb), hence 

n 

Ab{A){x) = B{A){{A-Lb){x)) = B{A){A{x)-Bx) = B{A){0) = ^Bfe_ii*"-'=(^). □ 
In addition to the notation given in relation (U), we also use 5k {B; t/) = [ 5^ (S; if) ... (i?; w) . 
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Lemma 3.2 Using the above mentioned notation and assuming [vi ... Vn]'^ ^K'^''^^ , we have: 

n n 

^l{B]^hijVj, . . . ^YynjVj) + 5l^^{B;vi, . . . ,Vn) = 5 k{B) ■ . (6) 
i=i 3=1 

Remark 3.3 The previous result can be described by: 

S^k{B-Bv)+5l+^{B;v) ^5k{B)-v,, l<z<n, (7) 
or simply by the following vector equation: 

5k [B- Bv) + Sk+i {B;v)= 5k [B) ■ v. (8) 

Proof of Lemma 13.21 Let e, £K"^^^ denote the column whose only nonzero entry is 1 in s*'' position. We 
also write Bj^s for s*'' column of the matrix B and [B]s for a square matrix of order n—1 obtained from B by 
deleting its s*'' column and row. According to the notation used in (j?]), let [B'^{Bis)]s stand for the matrix 
of order n—1 obtained from B by substituting s*'* column B^s in place of i*'' column, and then by deleting 
s*'' column and s*'* row of the new matrix. By applying linearity of Sk{B; ifj with respect to v, we have: 

n n 

5l{B;Bv) + 5l+,{B;v) = 5l{B;J2vsB^,)+5l^,{B;J2vses) 

71 11 n 

= Y^VsSk {B- B^,) +J2vs5l+, {B; e.) = {51 {B; B^^) +51+, {B; )) 

5=1 S— 1 S — 1 

n 

= v,{5l{B;B^,) + 5l+,{B;e,)) +J2vsi5k{B;B:,s) +5l+,{B;e.)). 

s=l 

First, we compute v,{51{B; B^,) +51+,{B-, e,)) ^v,{51{B) +5k{[B]:^) ^v^5k{B) . 

n 

Then, it remains to show that 'y^ys[5l.{B; B^s) +5l._^_i(^B; e^)) =0. 

s=l 

It suffices to prove that each term in the sum is zero, i.e. -64.^) (S; e^) =0, for s 7^ i. 

Suppose s ^ i. We now consider minors in the sum 5^5; B^^)- All of them containing s*'' column are equal 
to zero, so we deduce 

51{B;B^,) ^51{B\B^,)) ^5\{\B\B^,)],). 

If s ^ i, then each minor in the sum 5\^i{B;es) necessarily contains i*'* column and s*'* row. By inter- 
changing i*'* and s*'' column we multiply each minor by —1. We now proceed by expanding these minors 
along s*'' column to obtain —1 times the corresponding fc*'' order principal minors of matrix B'^{Bis) which 
do not include s*'* column. Hence e^) = —5\(\^B^(Bis)]s) and the proof is complete. □ 

In the following lemma we give some interesting correspondence between coefficients Bk of the matrix 
polynomial B{\) = adj(A/— _B) and sums of principal minors 5k+i{B;ifj , 0<k<n. 
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Lemma 3.4 Given an arbitrary column [vi ... Vn]'^&K^' 

Sk+i{B;vi,...,v„ 



it holds: 



Vl 

















= {-!)'' 6 k+i{B;v). 



(9) 



Sl+i{B;vi,...,i 

Proof. The proof proceeds by induction on k. It is being obvious for fc = 0. Assume, as induction hypothesis 
(IH), that the statement is true for fc — 1. Multiplying the right side of the equation Bk = Bk-iB + dkl, by 
the vector v, we obtain: 

BkV^Bk-iiBv) + dkV = {-l)''-'Sk{B;Bv)+dkV 
(IH) 

= i-l)'^-\6,{B;Bv)-5,v) =i-l)'^6k+i{B;v). □ 



Remark 3.5 The preceding Lemma [3.4l seems to have an independent application. Taking v — ej,l<j<n, 
we prove formulae (8) — (10) given by [1]. 

Theorem 3.6 (The main theorem for totally reduced systems - vector form) Linear system of first order 
A-operator equations ^ can be reduced to the system of n*'^ order A-operator equations: 



{AB{A))ix) = J2i-lf-'MB;A-'^{^)). 



(10) 



fc=i 



Proof. It is an immediate consequence of Theorem 13.11 and Lemma 13.41 

n n 

Ab{A){x) = Y,B,-i{A-\0)) ^Y.'^-l)'^'%{B;A"-'^{0)).n 

We can now rephrase the previous Theorem as follows. 

Theorem 3.7 ( The main theorem for totally reduced systems) Linear system of first order A-operator 
equations ([T]) implies the system, which consists of n^^ order A-operator equations: 



(As(A))(xi) = ^(-l)'^-Ml(i3;A"-n^)) 
fc=i 



(AB(A))(x,)-^(-l)'="MUi?;i'"-n^)) 



> . 



k=l 



(As(A))(x„) = ^(-i)'^-i(5:(S;A"-'=(^)) 



(11) 



fe=i 

Corollary 3.8 If we take A — 0, Theorem \'6 .7\ is equivalent to Cramer's rule. 

Remark 3.9 System (fTT|) has separated variables and it is called totally reduced. This system consists 
of n*'' order linear A-operator equations which differ only in variables and in the non-homogeneous terms. 

Remark 3.10 Formulae of total and partial reduction [7] could be of interest in the theory of linear control 
systems [2] and also in problems concerning systems invariants under the action of the general linear group 
GL{n,K), see [4] and [5]. 



Applications 

1° Assume that A is the differential operator on the vector space of real functions. Let us consider 
system ([1]) with initial conditions Xi{tQ) ~ Ci, 1 <i <n. The Cauchy problem for system ([T]) has a unique 
solution. Using form ([2]) , we obtain additional n — l initial conditions of i*'' equation in system ([TT1) : 



k=0 



l'=(^)(to)],, l<j<n-l, 



(12) 



where [ ]i denotes ■i*'' coordinate. Then each equation in system ([TT|) has a unique solution under given 
conditions and additional conditions (|12|) . and these solutions form a unique solution to system ([T]). Therefore, 
formulae ([TT]) can be used for solving systems of differential equations. 

Let us emphasize that we can extend this consideration to systems of difference equations and whenever 
the above method can be generalized. 

2° Now suppose that V is the vector space of meromorphic functions over the complex field and that A is 
a differential operator, i.e. A{'w) = -j^iw)- Let us consider system ([T]) under these assumptions. Recall that a 
function w &V \s differentially transcendental if it does not satisfy any algebraic differential equation. 
We prove the following theorem. 

Theorem 3.11 Let ipj he the only differentially transcendental component of the free column (p. Then for 
any solution x of system ([l} , corresponding entry Xj is also differentially transcendental function. 

Proof. We can see that the sum X]fc=i ^"'~'^('^)) is a differentially transcendental function. 

According to Theorem 2.8 proved by [6], the following equation: 



k=l 

yields that Xj IS 3( differentially transcendental function too. □ 
4. Examples 

Finally, we can illustrate our results with two special cases of n = 2 and ti = 3. 
Example 4.1 Consider the linear system of the first order A-operator equations, in unknowns xi,X2'- 

A{xi) = biixi + bi2X2 + ipi 

A{X2) = b2lXi + b22X2 + IP2 

Applying Theorem 13 . 71 the totally reduced system of system (fT3|) has the form: 

Ab{A){xi) = A^{xi) - SiAixi) + 62lixi) ^ 

bi2 



(13) 



A{ip2) b 



22 



bi2 

•^2 622 



Ab(A)(x2) = A^(X2) - bxA{x2) + ^2l{:X2) 

biiA{i^i) 

621 A{ip2 



bii 




b2i 


^2 



where we take into consideration only framed minors. 
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Example 4.2 Consider the linear system of the first order ^-operator equations, in unknowns xi,X2,X3: 

A{xi) = 611X1 + bi2X2 + bi3X3 + ipi 

A{X2) ^ b2lXi + b22X2 + b23X3 + ip2 /• (14) 

^ A{x3) = bsixi + 532X2 + bssxs + ips ^ 

Applying Theorem 13.71 the totally reduced system of system has the form: 

' Ab{A){xi) = A'ixi) - 6iA^{xi) + 62A{xi) - Sslix,) 

A^{ipi) bi2 bi3 



A'^{ip2) &22 623 
^^(^3) ^32 ^33 



^(■^i) bi2 


bi3 


A{(p2) b22 


b23 


A{(P3) b32 


&33 



As(A)(x2) 



A^X2) ~ SlA^(x2) + 52A{X2) ~ 53/(2:2) 

611 A'^iipi) bi3 



b2i 



631 A^iifs) 633 



fell 




bi3 


^21 


A{ip2) 


623 


&31 


Ai^3) 


^33 





bi2 


bi3 




b22 


623 


f3 


b32 


&33 




bii 


VI 


bi3 


b2i 


^2 


623 


b3i 


fa 


633 



As(A)(x3) = A^Xi) - SiA^ix3) + S2A{X3) - 63I{X3) 

bn bi2A^{ipi) 

&2I ^22 A'^{ip2) 
b31 632 A'^{(p3) 



fell 


fel2 A{ipi) 


fe21 


b22 A{(p2) 


fe31 


b32 A{ip3) 



fell 


fel2 


fl 


fe21 


fe22 


^2 


fe31 


fe32 


f3 



where we calculate only sums of the framed minors. 
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